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ABSTRACT 

We determine the initial mass function (IMF) of the "thin disk" by means of a direct 
comparison between synthetic stellar samples (for different matching choices of IMF, 
star formation rate SFR and depletion) and a complete (volume-limited) sample of 
single stars near the galactic plane {\z\ < 25 pc), selected from the Hipparcos catalogue 
(d < 100 pc, My < +4.0). Our synthetic samples are computed from first principles: 
stars are created with a random distribution of mass M* and age t* which follow a 
given (genuine) IMF and SFR(£* ) . They are then placed in the HR diagram by means 
of a grid of empirically well-tested evolution tracks. The quality of the match (synthetic 
versus observed sample) is assessed by means of star counts in specific regions in the 
HR diagram. 7 regions are located along the MS (main sequence, mass sensitive) , while 
4 regions represent different evolved (age-sensitive) stages of the stars. 

We find a bent slope of the IMF (using the Scalo notation, i.e., power law on 
a logarithmic mass scale), with T\ = —1.70 ± 0.1 (for 1.1A/ Q < M* < 1.6M Q ) and 
T2 = —2.1 ±0.1 (for « 4M > A/* > 1.6M Q ). In addition, comparison of the observed 
MS star counts with those of synthetic samples with a different prescription of the MS 
core overshooting reveals sensitively that the right overshoot onset is at M* — 1.50Mq. 

The counts of evolved stars, in particular, give valuable evidence of the history of 
the "thin disk" (apparent) star formation and lift the ambiguities in models restricted 
to MS star counts. A match of the evolved star counts yields a stellar depletion in 
the sample volume which increases with age (i.e., by apparent SFR(age) / SFRo). 
A very good match is achieved with a simplistic diffusion approximation, with an 
age-independent diffusion time-scale of r<jif = 6.3 • 10 9 yrs and a (local) SFR = 
2.0 ±0.1 stars (with > 0.9M Q ) formed per 1000 years and (kpc) 3 . We also discuss 
this "thin disk" depletion in terms of a geometrical dilution of the expanding stellar 
"gas", with H z (t„) oc aw(t*)- This model applies to all stars old enough to have 
reached thermalization, i.e. for f» > 7 • 10 8 yrs and H z > 230 pc. It yields a column- 
integrated (non-local) "thin-disk" SFR co i which has not changed much over time (< 
30%), SFR col w 0.82 * pc^Gyr" 1 (i.e., stars with M» > O.9M ). 

Key words: Stars: evolution - Stars: late-type - Stars: mass function - Galaxy: 
kinematics and dynamics - Galaxy: disk - Solar neighbourhood 



1 INTRODUCTION 

The case of the IMF has spurred a vast amount of literature 
over the past 3 decades. We would like to mention only the 
classic work and reviews by Miller & Scalo (1979) and Scalo 
(1986), together with Scalo (1998) and Kroupa (2001) for 
reference to the more recent work in the field. The classical 
approach to derive the IMF of field stars in the galactic disk 
follows a line of steps: (i) derive the stellar distribution in ab- 
solute magnitude from the observed distribution in apparent 
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magnitude, requiring the respective vertical space-density 
scale-height, and convert that into the luminosity function 
LF for MS stars, (ii) apply a theoretical mass- luminosity 
(m-L) relationship for obtaining the present-day mass func- 
tion (PDMF), and (iii) model the IMF on the PDMF, for 
a given SFR history. As discussed by Scalo (1998), each of 
these steps introduces specific uncertainties so that the case 
of the IMF is far from being closed. Scalo (1998) and Kroupa 
(2001) even argue that there are systematic variations in the 
IMF, especially between the field star IMF and the various 
IMF's derived from galactic stellar clusters. 

The availability of good parallaxes, as required for a 
precise field star LF, is mutually exclusive with the neces- 
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sity of a large stellar sample. Therefore, early work on the 
LF had to rely on a statistical approach, based on the appar- 
ent magnitudes, or on photometric parallaxes (e.g., Gilmore 
& Reid 1983). These approaches link the question of the 
IMF to that of the spatial stellar distribution, i.e., perpen- 
dicular to the galactic plane. Most notably, Gilmore & Reid 
(1983) found that (i) the stellar spatial density scale heights 
increase with lower masses, and that (ii) the global spatial 
distribution could be described by a (denser) Galactic "thin 
disk" plus an extended Galactic "thick disk" . Especially the 
first point is of importance when deriving a field star IMF, 
as we will show below. 

The great interest in the IMF is fuelled not only by the 
attempts to understand and model the star formation pro- 
cess (es). Finding the IMF is also closely linked to the ques- 
tion of galactic star formation history, i.e., whether and how 
the SFR has varied over time. Both IMF and SFR(t) are im- 
portant ingredients of galaxy models since these quantities 
determine such important aspects as the chemical evolution 
and the luminosity evolution of a galaxy (see, e.g., Dwek 
1998, Gilmore 1999, Hensler 1999, Pagel 2001a,b). Further- 
more, the IMF of more massive stars is of particular interest 
for the formation of star-burst regions and stellar clusters, 
and for the evolution of young galaxies. Our Galaxy would 
make an ideal test case for detailed galaxy evolution models, 
if we could determine the IMF and SFR(t) with sufficient ac- 
curacy. Already, we have quite a good account of the stellar 
population(s), the (local) kinematics and chemical composi- 
tion as a function of age. 

A general problem in this context is the inherent am- 
biguity of any interpretation that is restricted to MS stars: 
Towards lower stellar masses, the mean age of a MS star 
increases accordingly. Consequently, assuming a different 
time-dependence of the rate of star formation (i.e., increas- 
ing, or decreasing with age) has the same effect on the ex- 
pected numbers of stars on the lower MS as changing the 
IMF accordingly. This degeneracy of SFR(t) and lower IMF 
(see, e.g., Binney et al. 2000) allows a variety of false com- 
binations of IMF and SFR(t) to match the PDMF of MS 
stars. It is therefore important to lift this degeneracy by 
extracting additional, specifically age-sensitive information 
from the HR diagram which is buried in the properties and 
numbers of the evolved stars. 

Matters are further complicated by the galactic disk 
dynamics. Relevant properties of the star formation process 
are the genuine IMF and SFR at the time of the star for- 
mation and in a defined volume on the galactic plane - i.e., 
within the vertical extent of the star formation sites. How- 
ever, present-day counts of stars in a solar neighbourhood 
volume translate into an apparent SFR which decreases to- 
wards larger age (Schroder & Sedlmayr 2001). This deple- 
tion of stars in the galactic plane is caused by a dilution of 
the "heated stellar gas" which is expanding vertically (see 
below). If not accounted for, this effect would make a spuri- 
ous, systematic difference to, e.g., the steepness of any IMF 
derived from a volume-limited stellar sample. Consequently, 
the above-mentioned classical studies of IMF and SFR have 
used the invariant column-integrals. Their dependence on 
stellar mass and age is equivalent to that of the genuine 
IMF and the depletion-corrected SFR in the volume. This 
column-integral approach, however, has the disadvantage of 
being no more accurate than the required and presently 



quite limited knowledge of space-density scale-heights as a 
function of stellar age. 

For this reason, we use the direct approach to star for- 
mation, based on the stellar records in a well-defined (local) 
volume, and an interpretation in terms of a matching syn- 
thetic sample. Synthetic stars are randomly placed in a well- 
tested evolution grid, and their numbers are (statistically) 
defined by the given (genuine) IMF and (apparent) SFR ap in 
the volume. In previous studies we have shown that counts 
of specific synthetic stars compare well with the respective, 
observed present-day star counts (Schroder 1998; Schroder 
& Sedlmayr 2001). The potential of this approach lies in its 
simplicity. We do not need to determine a LF, nor to use a 
(m-L) relation - this is automatically accounted for by our 
synthetic sample creation, as is the mass-dependent fraction 
of post-AGB stars which are left out by a MS-based PDMF. 
The genuine SFR(t) in the volume is then derived from a 
depletion description (see below), for which a maximum of 
age-sensitive information is derived from the stellar sample. 

Studies based on a similar approach have been under- 
taken by Bertelli & Nasi (2001), who, in addition, account 
for a variety of metallicity (see also Bertelli et al. 1992), and 
by Binney et al. (2000). However, these authors did not ex- 
ploit the specific, age-sensitive information provided by the 
counts of evolved stars. In fact, the models of Bertelli & 
Nasi (2001) show an excess of the number of evolved stars 
by (relative to the MS stars of same initial mass) a factor 
of 1.5, while the problem of stellar diffusion away from the 
galactic plane has not been addressed. 

Consequently, the work presented here pays special at- 
tention to three points which appear crucial for an im- 
provement over previous work: (1) In section 2, we describe 
a large, volume-limited sample of stars from the Hippar- 
cos catalogue, which we selected for good proximity to the 
galactic plane (\z\ < 25 pc) in order to avoid any selective 
under-representation of young, massive stars at larger \z\. 
(2) In section 3, we study the stellar depletion in the galac- 
tic plane. A simplistic but well-matching diffusion approxi- 
mation with an age-independent diffusion time is compared 
with a more physical model. The latter describes the de- 
pletion in terms of geometrical dilution caused by the stel- 
lar kinematics: stars old enough (£* > 7 ■ 10 yrs) to have 
"thermalized" are expanding to larger scale heights exactly 
as their vertical velocity distribution increases with age. We 
find a semi-empirical prescription of the stellar space-density 
scale-heights in agreement with the available observational 
evidence. (3) In section 4, we describe the computation of 
our synthetic samples and how we assess a match with the 
observed stellar sample. Special consideration is given to the 
counts of evolved stars, in order to constrain the right choice 
of IMF, SFR(t) and depletion as well as possible. 



2 COMPLETE SAMPLES OF SINGLE STARS 
ON AND NEAR THE GALACTIC PLANE 

In principle, we use the same volume- limited (i.e., d < 
100 pc) sample of single stars, based on the Hipparcos cat- 
alogue, as for earlier work (Schroder 1998), but solving the 
problem with consistency: In comparison to the very local 
sample with d < 50 pc (used in Schroder & Sedmayr, 2001), 
the lOOpc-volume provides a comfortably larger statistical 
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Figure 1. The observed HR, diagram of all single stars which are within 25 pc to the galactic plane, less than 100 pc away, and brighter 
than My = 4.0, as obtained from the Hipparcos catalogue entries. Also shown are the regions used for characteristic star counts (see 
text), 7 along the main sequence and 4 for evolved stages. 



stellar basis, but it also comprises larger distances z from 
the galactic plane which already make a difference to the 
PDMF of massive stars. Therefore, this investigation now 
discriminates against all stars with \z\ > 25pc within 100 pc, 
which still leaves us with larger, statistically more signif- 
icant star counts than from the 50pc-sample. Star counts 
from neighbouring space are used, where possible, to derive 
space-density scale-heights for specific groups of stars (see 
Table 1). 

For computing the z values, we adopted an offset of 
the Sun of 15pc north of the galactic plane (Cohen 1995), 
which intersects the spherical sample volume around the Sun 
(R = lOOpc, V = 4.19 ■ 10 6 pc 3 ) off-centre. In order to keep 
the star counts of the neighbouring space reasonably large, 
we devide the spherical volume into slabs of (A) all stars 
within \z\ < 25pc (Va = 1.50 ■ 10 6 pc 3 ), (B) all stars between 
25pc< \z\ < 50pc (V B = 1.31 • 10 6 pc 3 ), and (C) those stars 
in the polar sections of the lOOpc sphere, up to \z\ < 85pc 
and 115pc, respectively (Vc = 1-38 • 10 6 pc 3 ). 

Stellar population densities in the HR diagram are as- 
sessed by means of star counts in 7 regions along the main 
sequence (MSI to MS7, sensitive to the initial mass function 
IMF), and in 4 regions representative of different evolution- 
ary stages and stellar ages, containing information on the 
star formation history SFR(t„) and the diffusion losses over 
time: the Herzsprung gap (HG), the K giant clump (KGC), 



the lower giant branch (LGB) and stars with cool winds 
(CW), with border-lines as shown in Fig. 1. The observed 
HG counts are particularly affected by any unrecognized bi- 
naries with a composite colour and must be considered as 
upper limits, only. 

Since the Hipparcos catalogue is complete at least for 
stars brighter than V = 7.3 (in some regions of the sky even 
down to V = 9, see Perryman et al. 1997), completeness of 
volume-limited samples is automatically given for all stars 
within 50pc which are brighter than My = 4. Towards lOOpc 
distance, however, star counts of regions MS6 and MS7 (with 
My = 3.0 to 4.0) would be somewhat incomplete. By apply- 
ing the (unaffected) star count ratios in the 50pc sample (as 
from Schroder & Sedlmayr 2001), especially to MS6, MS7 
and the total number of stars (TN) as a multiple of MS5, we 
can make a good estimate of the missing fractions. For the 
stars within \z\ < 25pc of the galactic plane, we find correc- 
tions for the counts of MS6 by 4%, of MS7 by 25%, and of 
TN by 8%. These figures need to be added for a comparison 
with the synthetic samples discussed in section 4. 

Normalizing the characteristic star counts to their re- 
spective volumes (A, B and C as given above), we derive 
approximate, local scale heights, as listed in the bottom line 
of Table 1. In order to handle any intrinsic fluctuations of the 
star counts, we applied a comparison of the ratios AT*,a/W»,b 
and Ar»,A/W*,c with volume integrals over exponential den- 
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Table 1. Characteristic star counts in the space slabs A, B and C, for different regions in the HR diagram along the main sequence 
(MS) and for different groups of evolved stars (see Fig. 1 and text). The bottom line gives the approximate scale heights of the respective 
stellar number densities. 
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HG 
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LGB 


CW 
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[0-25] 


38 


106 


317 


327 


484 


712 


850 


3348 


< 49 


205 


157 


16 


B 


[25-50] 


13 


79 


227 


261 


365 


581 


743 


2736 


< 53 


171 


146 


21 


C 


[50-115] 


20 


63 


212 


260 


374 


644 


835 


2921 


< 69 


175 


156 


22 




Hz,* ~ 


60: 


135 


160 


320: 


270: 


> 400 




> 400 




> 400 







sity distributions for different scale heights. For the sake of 
a more realistic description, we assume a constant space- 
density around the galactic plane for |z| < 10 pc. 

The small space-density scale-heights of the stars in 
groups MSI to MS3 confirm that the young and more mas- 
sive stars are much more concentrated to the galactic plane, 
and that their correct assessment requires a volume restric- 
tion in z. For MS6, MS7 and LGB, in fact, stellar number 
densities increase slightly from B to C. This local deficiency 
of stars appears to correspond with the low-z-structure 
found by Kuiken & Gilmore (1989). Hence, the // z -values 
stated in Table 1 should be larger than what is typical for the 
galactic thin disk over a larger z-range. Clearly, a complete 
sample covering a larger volume, as will be available from 
the future GAIA (e.g., Perryman 2002) and DIVA (Bastian 
et al. 1996, Seifert et al. 1998) missions, is highly desirable. 



3 VERTICAL SPACE DENSITY STRUCTURE 
AND "THIN DISK" DYNAMICS 

A lot of observational work has been done on the stellar 
space-density scale-heights H z vertical to the galactic plane. 
Galactic structure studies on a wide range of stars have come 
up with the concept of a "thin" disk and an old "thick" 
disk, with scale heights of « 250 pc to 400 pc (thin disk, 
e.g., Gilmore & Reid 1983, Kent et al. 1991, Ojha et al. 
1996, Siebert et al. 2002) versus 700 pc to 800 pc (thick 
disk, e.g., Ojha et al. 1999, Robin et al. 1996, Kerber et al. 
2001, Phleps et al. 2000). However, a single scale height does 
not describe the "thin disk" stellar distribution adequately. 
Rather, each particular class of stars has a different scale 
height which increases with the respective mean stellar age. 
Evidence comes from selective star counts which focus on 
narrow age ranges, though specific values may still be quite 
uncertain: 

(i) From an interstellar extinction study, Vergely et 
al. (1998) find a scale-height of the local ISM density and 
clouds, the regions of star formation, to be between 35 pc 
and 55 pc. 

(ii) For young and massive stars, scale-heights between 
25 pc and 65 pc have been found (Reed 2000, for OB stars) , 
as well as 45 pc (Conti & Vacca 1990, WR stars). 

(iii) From a study of all F-type stars within a distance of 
50 pc, which (on average) are half-way through their main- 
sequence life, Marsakov & Shelevev (1995) find a scale-height 
of 160 pc. 

(iv) For thin disk K dwarfs, Kuiken & Gilmore (1989) 
find a scale-height of 249 pc for z > 300 pc. Nearer to the 
galactic plane, their density data show pronounced low-z- 



structure. A rough average would yield a significantly larger 
scale- height. 

(v) For the galactic disk subgiants, Mendez & van Al- 
tena (1996) find a red-giant scale-height of 250 pc. Many 
of these are He-burning K-giant-clump stars, represented by 
our group KGC. 

(vi) Jura & Kleinmann (1992a, b) studied K and M- 
type (O-rich) variable AGB stars and Miras. For periods of 
P — 300. ..400 days, indicative of a larger luminosity and, 
hence, less small mass (around 1.1, possibly 1.2 M Q ) and 
lesser age, they find a typical scale-height of 250 pc. 

(vii) For the supposedly less luminous, and therefore 
presumably slightly older group of stars with P — 200. ..300 
days, the same authors estimate 500 pc. In addition, Jura 
(1994) finds a scale-height of 600 pc for the O-rich, short- 
period Miras, which he estimates at a mass as small as about 
1.0M Q . 

(viii) for RR Lyrae stars with a moderate metal under- 
abundance ([Fe/H]>-1) Layden & Andrew (1995) find a 
scale height of 700 pc. As expected from their abundance, 
this places them at the transition from thin to thick disk. 

A similar range of scale-heights emerges from a study 
of LMC foreground stars, grouped by different spectral class 
and metallicity, by Oestreicher & Schmidt-Kaler (1995). 
These authors find 105 pc for early-type stars, to over 500 pc 
for the metal-poor, late-type giants. 

In line with these observations, an analytical description 
of the stellar dynamics suggests age-dependent scale heights 
which increase in proportion with the vertical stellar velocity 
distribution a w : 

After a few oscillations at right angles to the Galactic 
plane, stars reach a state of dynamical equilibrium governed 
by the potential gradient at right angles to the plane and 
by the velocity dispersion (Oort 1932). With the progress 
of time, the velocity dispersion increases from the effect of 
encounters with interstellar clouds; alternatively, the oldest 
stars, belonging to the Thick Disk, are born with a high 
velocity dispersion, which they still retain, either as a re- 
sult of being born from collapsing gas clouds that had not 
yet undergone complete dissipation of their potential en- 
ergy (Eggen, Lynden-Bell & Sandage 1962; Burkert, Truran 
& Hensler 1992), or because of merger events in the early 
history of the disk (Robin et al. 1996; Pilyugin 1996). Either 
way, the velocity dispersion increases with age (Wielen 1977; 
Lacey 1984; Sommer-Larsen & Antonuccio-Delogu 1993) re- 
sulting in an increase in the width of the density distribution 
as a function of z. 

By Jeans's theorem, the distribution function in phase 
space in a steady state is a function of isolating integrals of 
the motion (Jeans 1915, Lynden-Bell 1962). In the approx- 
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imation where z-motions are considered as uncoupled from 
motion in the plane, the only isolating integral is the vertical 
energy 



1 2 

2 W +' 



(1) 



and the distribution function (assumed to be gaussian 
for an isothermal stellar population) is 



f{z,w) = v(z)- 



-(w 2 +2</>)/2o 



'WW (2) 

where v is the volume density at height z, w the ve- 
locity along the z-axis and a w the corresponding velocity 
dispersion. Consequently the total volume density is 



u( Z ) = r 

J — c 



f(z, w)dw — i/(0)e 



-*/< 



(3) 



Kuijken & Gilmore (1989) have given an expression for 



4>{z): 



<f> = 27tG[E (V z 2 + D* - D) + PcS z% 



(4) 



where E is the mass surface density of the disk, D a 
length scale defining the transition from the linear regime 
K z oc z to K z — const, and p e a is the contribution of the 
dark halo to the local density of gravitating matter. This 
expression can be expanded (for z 2 < D 2 ) as a binomial 
series 



by Siebert, Bienayme and Soubiran (2002), among others. 
Eo/D is rather well determined at 0.18 Mq pc -3 and p e ff ~ 
0.02 is a relatively small correction, but there is a degeneracy 
in the determination of Eo and D. Siebert et al. favour D w 
400 pc, by contrast to Kuijken & Gilmore's value of only 180 
pc. The difference is crucial to the range of validity of Eqs 6, 
8 and 9. Taking D = 500 pc, the quadratic approximation 
to <f> is accurate to better than 10 per cent up to z = 350 
pc, which would make it applicable to virtually all thin-disk 
stars in our neighbourhood, whereas with D = 180 pc, the 
situation would be very different. The actual scale heights 
of the more numerous, less massive stars presented above, 
however, suggest that the approximation is indeed valid for 
our sample. 

In the approximation of the density distribution viz) 
being an exponential function of z with a scale height H z , 
we have the simple relation N = 2- H z ■ v(0). In combination 
with Eq (9) it then follows that H z (t*) expands oc a w (t*) 
and u(0) dilutes oc l/a w (t*) with age and increasing vertical 
velocity dispersion. 

For a w (t*) as a function of stellar age, we take the 
simple (constant-diffusion) solution of Wielen (1977), a w — 
(a w + b w ■ t,) - 5 . For U > 7 ■ 10 8 yrs, we find that Wielen's 
empirical data points, as well as those of Sommer-Larsen & 
Antonuccio-Delogu (1993), can be well enough represented 
by a w = 72 (km/s) 2 , b w = 72 (km/s) 2 /10 9 yrs. 

Finally, with H z oc o w (t*) and in consideration of the 
above listed empirical stellar space density scale heights, we 
have adopted the following prescription (in pc) as a function 
of age (in yrs) : 



2ttG 



= 2ttG 



% . . , 1 z 2 1 z 4 
So ' D| 2^" 8D* + 



+ PeffZ 



/Eo_ 
\2D 



+ PeS Z 



(5) 



The first term, <j> oc z 2 , corresponds to simple harmonic 
motion where from Eq 3, 



Kz) 

1/(0) 



where 



A 2 = >,(; 1 Y 



(^+Peff) 



(6) 



(7) 



and in this approximation the column density of a tracer 
population by number is related to the central density by 



jV 



v(z)dz = v{0)-j£^/ir, 



(8) 



177- (1.0 + ^/10 9 ) ' 5 , t*>7-10 8 



(10) 



In particular, scale heights range from around 265 pc 
(about the average half-thickness of the thin disk according 
to Ojha et al. 1996, reached around an age of t* w 1.2 • 
10 9 yrs) to about 560 pc for the very oldest (9 • 10 9 yrs) 
thin-disk stars. 

The present-day, young stars are formed very near the 
galactic plane, with the scale-height of the ISM clouds 
{H z fi w 45 pc, see above), and drift away with increas- 
ing age. As long as their kinematics are not thermalized 
(for t* < 7 ■ 10 s , H z < 230 pc), we prescribe their scale- 
heights mainly by a quadratic superposition of the initial 
scale-height H Zy o and a term which is linear in age (in years), 
as: 



H z ,, = (45 2 + (t„/3.105 • 10 6 yrs) 2 ) ' 5 , U < 7 ■ 10 s 



(11) 



The resulting values are in reasonable agreement with 
the scale-heights obtained directly from the corresponding 
age-groups of Hipparcos stars (i.e., MSI, MS2, and MS3 in 
Table 1). 



1/(0) 



A 



-N, 



(9) 



i.e. we have a gaussian distribution with an effective 
scale height proportional to o w ; in a higher approximation 
to <f>(z), Eq (3) must be integrated over z numerically. 

The values of the parameters Eo, D and p e ff have been 
discussed by Kuijken & Gilmore (1989) and very recently 



4 "THIN DISK" DEPLETION BY DILUTION: 
EXPANSION OF THE STELLAR "GAS" 

Near the galactic plane, the vertical gravitational net force is 
small enough to be neglected. In an earlier study (Schroder 
& Sedlmayr 2001), we therefore suggested a very simplistic 
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model in which depletion was prescribed as a diffusion ef- 
fect, with a single diffusion-time scale Tdifi for all stars (now 
referred to as model a). A value of 6.3- 10 9 yrs, which we here 
confirm to give a good match of the observed star counts, 
leads to a 50% depletion of stars aged 4 • 10 9 yrs. To give an- 
other example of the significance of this depletion: the ratio 
of less massive giant stars (aged around 5 • 10 9 yrs) over their 
MS counterparts (of average age 2.4 • 10 9 yrs) is reduced by 
a factor of 1.5 (which is of the order of the mismatch found 
with synthetic samples by Bertelli & Nasi, 2001). 

For the IMF, we use the conventional power-law nota- 
tion (Scalo, 1986) that counts on a logarithmic mass scale: 

IMF(M„) = dn«/dlogM« oc Ml (12) 

To allow for a bent slope, we use a Ti for masses smaller than 
I.GMq, and a F2 for larger masses. The mass-range of this 
study is limited by the observed sample to about (effectively) 
AMq at the upper, and to about I.IMq at the lower end. 

Our principal star formation observable is the appar- 
ent SFR ap (i*) in the sample volume. As mentioned before, 
it appears to have been much lower in the past than at 
present (SFRo) due to a dilution effect of the expanding 
stellar "gas". Therefore, we introduce a depletion factor of 
fde P i(i*) = SFR ap (i*)/SFR(i») relative to the genuine star 
formation rate SFR(t») in the volume, at the time of the star 
formation. In addition, we consider possible slow changes in 
the genuine SFR by an exponential factor, relative to the 
present-day star formation rate SFR , defined on the re- 
versed time-scale of the stellar age t* (in years): 

SFR ap (^) = Fdepi^*) ■ SFRo • e ^«/ 1()9 (13) 

Our values of the SFR (in the volume) are given in stars 
formed per year which have a mass larger than 0.9Mq , in the 
local sample volume A (see previous section) of 1.5 • 10 6 pc 3 . 
For our simplistic diffusion model (a), the depletion factor 
fdcpi (for any given number n* of similar stars of a specific 
age and mass) is itself a f*-dependent, exponential term (in 
Schroder & Sedlmayr 2001, both these terms were simply 
merged into one). 

A more physical description of the depletion is given in 
terms of a geometrical dilution of stellar space density into 
the column (model b). According to Eq 9, the depletion fac- 
tor is then proportional to the vertical velocity distribution 
or, equally, to the scale-height of a group of stars of a given 
age: 

F dep i(U)=H ie{ /H z (U) (14) 

This model (b) applies, once (1) the stellar "gas" has 
thermalized (i.e., if t* > 5 • 10 8 yrs or more), and (2) if the 
depletion by dilution into the column is not compensated 
for by any fill-in from radial mixing in the disk. However, 
for the past 7 • 10 yrs, significant fill-in is evident: From 
this more recent time-span we would expect a dilution by a 
factor of about 5 (from the scale- height of the star- formation 
layer, about 45 pc, to H Z (U = 7 - 10 8 ) = 230 pc. Instead, the 
respective Fd ep \(t t ) required by a matching synthetic sample 
(see below) is about 0.85 or even closer to 1. 

To match this depletion at the onset of geometrical di- 
lution (7- 10 8 yrs), H TC { must be 195 pc. For all younger (i* < 



7-10 yrs) stars we use the simple exponential approximation 
of the depletion factor this time with Tdu = 4.5 • 10 9 yrs), 
guided by a smooth transition with the depletion derived 
from dilution. 

For a quite limited volume around the Sun, a signif- 
icant fill-in over the more recent past must be expected: 
At present, it is located right in between two spiral arms. 
Consequently, the present-day, local SFR is much below the 
global, average SFR for the "thin disk" . As we look towards 
older stars, radial mixing changes the observed SFR from 
a local (for the youngest stars) to a non-local (after 2 to 3 
orbital periods) quantity. 

A quantity more relevant for the physics of the galactic 
disk is the column-integrated star formation rate, which re- 
lates to the apparent local SFR in the volume, our principal 
observable, by: 

SFR co i(i*) « 2-H z {U) -SFR ap (15) 

This column-integrated SFR is then related to the real 
local star formation in the volume (allowing here for a pos- 
sible exponential dependence on age) and the subsequent 
depletion, simply by: 

SFR co i(t«) w 2 • H Z {U) ■ F depi {U) • SFR • e 1 '*' 10 " (16) 

The SFR co i(t*) derived from our model first increases 
with look-back time (driven by the above-mentioned fill-in 
from radial mixing which causes scale-heights of the younger 
stars to increase steeply with age), until from an age of 
7- 10 8 yrs it then settles on a disk-averaged value (see below). 
This is about 4 times larger than that of the low present-day, 
interarm star formation. The same effect would be observed, 
if the effective star formation layer had just been more ex- 
panded in the past (as noted by Binney et al. 2000). 



5 SYNTHETIC STELLAR SAMPLES AND 
THE IMF IN THE GALACTIC PLANE 

— The synthetic samples: Our synthetic stellar samples 
are computed from first principles: stars are created with a 
random choice of mass M* and age t* within a distribution 
which is prescribed by the choice of IMF for stellar masses, 
and SFR(i*) for the ages. These synthetic stars are then 
placed on a fine-meshed grid of individually computed evolu- 
tion tracks, in order to obtain their present-day distribution 
in the observational HR diagram (see also Schroder & Sedl- 
mayr 2001). We computed the evolution tracks with a fast 
but well-tested code, developed by Peter Eggleton (Eggle- 
ton 1971, 1972, 1973). This code has been updated with 
up-to-date opacities and colour tables. Its performance and 
the right choice of mixing length and overshoot prescription 
has been tested very critically in several studies. Respective 
evolution models have been compared with a number of real 
MS and giant stars in well-observed eclipsing binaries, and 
isochrones have been compared with stellar cluster data (see 
Schroder et al. 1997, and Pols et al. 1997, 1998). 

Because of the limited observational data, we did not 
explicitly account for the effects of a moderate range in 
chemical composition, which mostly contribute to a more 
smeared-out distribution of the stars in the synthetic HR 
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Table 2. Characteristic star counts Ni of different synthetic samples and their deviations from observed values, in units of the statistical 
variation y/n^obsi P ms the mean deviation |<5| av (double weight is given to KGC and LGB counts). Under the best-matching, simple 
diffusion model, synthetic samples based on different evolution grids and thin disk ages are listed, followed by those derived with depletion 
by (mainly) dilution. A constant SFR (7 = 0) gives a better match than SFRs decreasing from an initially 30% (7 = 0.03) and 70% 
(7 = 0.06) higher value. 



Model 


MSI 


MS2 


MS3 


MS4 


MS5 


MS6 


MS7 












HG 


KGC 


LGB 


CW 


|<5|av 


Obs. 


38 


106 


317 


327 


484 


740 


1068 












(< 49) 


205 


157 


16 




T dif = 6.29 


37 (-o.io) 


109 (+0.29) 


315(-°- 14 ) 


325 (-0.13) 


495 (+0.48) 


73g (-0.04) 


1077 (+0.27) 












44 


209 (+0.28) 


1 49 (-0.61) 


18 (+0.40) 


0.30 


ovo=1.8 


40<+0-25) 


108 (+o.i7) 


317 (+o.o) 


260(- 3 7 ) 


472 (-0.50) 


770 (+l.l) 


107 g(+0.23) 












44 


l 97 (-0.59) 


l 79 (+i-8) 


20 (+i.o) 


0.98 


ovo— 1 .6 


4Q(+0.2B) 


10 g(+0.27) 


onc( + 0.44) 


3n (-0.87) 


451 (-1.5) 


764 (+0.87) 


1062( _o ' 17 ) 












42 


206(+ - 11 ) 


162 (+0.36) 


ig(+0.65) 


0.50 


ovo=1.4 


37 (-o.i3) 


112 (+0.5S) 


322 (+0.29) 


318 (-0.50) 


539 ( + 2.5) 


75l(+o.4i) 


1064(-°- 13 ) 












39 


200 (-0.36) 


147(-o.77) 


22 (+1.5) 


0.69 


^9max 10 


38 (-0.06) 


107(+ - 14 ) 


3n (-0.33) 


330 (+0.19) 


507 (+1.05) 


716 (-0.89) 


1067 (-0.04) 












39 


l 99 (-039) 


159 (+0.18) 


18 (+0.44) 


0.36 




41 (+0.45) 


112 (+0.58) 


317 (-0.02) 


324 (-0.17) 


498 (+0.65) 


730 (-0.36) 


1079 (+0.34) 












48 


19l(-0.98) 


142 (-1.20) 


16 (+0.10) 


0.59 


Tdifoo, Y\ — 0.5 


42 (+0.68) 


10 g(+0.16) 


311 (-0.34) 


329 (+o.io) 


519 ( + 1.59) 


744 (+0.15) 


1073 (+0.16) 












56 


243 ( + 2.6) 


220 (+5.1) 


3 l(+3.7) 


1.9 


Dil.(t. > 7 10 8 ) 


39 (+o.io) 


114 (+0.74) 


331 (+0.78) 


323 (-0.21) 


500 ( + 0.72) 


725 (-0.57) 


1071 (+0.09) 










40 


206 (+0.21) 


164 (+0.59) 


21 (+1.35) 


0.49 


r 2 - 1.9 


45 (+i.i) 


116(+ 1 -°) 


322 (+0.29) 


3H(-0.87) 


49 l(+0.30) 


691 (-i.8) 


1061 (-0.22) 












40 


l 96 (-0.66) 


164 ( + 0.51) 


23 (+1.65) 


0.79 


r 2 - 2.0 


41 (+0.52) 


125 (+1.9) 


330 (+0.72) 


331 (+0.22) 


48 g(+0-25) 


700<- 1 ' 5 ' 


1087 (+0.58) 












39 


205 (+o.oo) 


162 (+0.43) 


21 (+1.20) 


0.64 


r 2 - 2.2 


34 (-0.58) 


115 (+0.89) 


322 (+0.30) 


319 (-0.46) 


503 (+0.87) 


727 (-0.48) 














45 


193 (-0.84) 


166 (+0.075) 


23 (+1.8) 


0.82 


r 2 - 2.3 


28 (-1.6) 


n2 (0.60) 


317 (+o.oi) 


312 (-0.84) 


505 (+0.97) 


764 (+0.90) 


H 24 (+1.7) 












43 


208 (+0.24) 


17l(+i-i) 


21 (+1.3) 


0.89 


Ti - 1.60 


39 (+o.i6) 


12 3 (+i.7) 


327 (+0.59) 


331 (+0.24) 


507(+ 1 -°) 


713 (-1.0) 


1107(+ 1 - 2 ) 












43 


204 (-0.04) 


l 59 (+o.i9) 


22 (+1.55) 


0.66 


Tj - 1.80 


41 (+0.4B) 


ll (+0-43) 


315 (-o.ii) 




47g (-0.23) 


727 (-0.46) 


1Q77 (+0.28) 












46 


195 (-0.73) 


l 66 (+0.72) 


23 (+1.7) 


0.64 


70.03 


3g (+0.03) 


113 (+0.70) 


326 (+0.52) 


316 (-0.61) 


47 g(-0.22) 


709 (- 1 1 ) 


m5 (+14) 












43 


206 (+0.08) 


l 75 (+i.4) 


20 (+0.91) 


0.72 


70.06 


38 (+0.03) 


ll (+0-37) 


318 (+0.04) 


326 (-0.08) 


47 g(-0.26) 


72l(-o.69) 


n23 (+1.69) 












46 


214 (+0.63) 


203 ( +3.6) 


25 ( + 2.3) 


1.16 



diagram. Instead, we used a random smearing prescription 
(see Schroder 1998) which would match the appearance of 
the observed HR diagram. It comprises the residual distance 
uncertainties in the Hipparcos parallaxes, photometric er- 
rors and some variation in metallicity. 

Table 2 lists the characteristic star counts of different 
synthetic stellar samples, starting with the best match which 
we obtained with the simple diffusion approximation (model 
type a), when using an IMF with Ti = —1.65 and T 2 = 
—2.05, an overshoot-onset (ovo) at 1.5OM0 (see below) and a 
constant (genuine) SFR (7 = 0.0) in the volume. Obviously, 
this simplistic model gives the best empirical description of 
the depletion. Therefore, we use it as a reference for all other 
synthetic samples. In all our models (with the exceptions 
of £9maxl0 and f9max8, see below), star formation starts 9 



billion years ago, which is a commonly adopted age of the 
thin galactic disk (Hernandez et al, 2001). 

— Evaluation of match quality: In order to decide the 
match quality of any synthetic sample on a quantitative 
scale, we compute the average (|<5| a v, listed by the last col- 
umn in Table 2) of all individual deviations from the ob- 
served star counts rii, in units of their individual statistical 
variations ^/ru. We prefer the choice of a simple (linear) aver- 
age over the rms error because not all individual deviations 
are independent of each other. Also, a double weight has 
been given to the deviations of the age-sensitive KGC and 
LGB counts for balancing these with those of the seven MS 
counts. Of the observed complete sample, counts for MS6 
and MS7 have been corrected for some incompleteness, see 
previous section. The synthetic star counts were obtained 
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from 5 times larger samples, to reduce their statistical con- 
tributions to the deviations. 

Besides their average, the individual 5-values are 
equally important for the choices of best-matching parame- 
ter. In fact, each of the parameters discussed below affects 
only certain counts in a very specific way, while the corre- 
sponding change in |i5| av appears to be much less significant. 
In this way, degeneracies with multiple parameter changes 
are avoided. In Table 2, we therefore highlighted individual 
mismatches by bold-face printing of the respective S- value. 

— Onset of overshooting on the MS: Because of the now 
larger (than in previous work) and therefore more critical 
star counts along the MS, we had first to reconsider the 
exact onset of MS overshooting by comparing the effects 
of different evolution grids on the MS star counts. For this 
comparison, we used the simple diffusion model because it 
would not introduce any artificial count-change along the 
MS (while the artifically abrupt onset of dilution at t* > 7 ■ 
10 8 yrs in our model type (b) does, marginally so). The best 
match (i.e., lowest mean deviation of its star counts from the 
observed record) is then achieved with an overshoot-onset 
(ovo) at I.50M Q , see Table 2. This is in good agreement with 
Pols et al. (1998) but slightly different from our previous 
choice of 1.8M Q . 

By contrast, the star counts of synthetic samples com- 
puted with ovo = 1.8, ovo = 1.6, and ovo = 1.4 (all pa- 
rameters for IMF and SFR left unchanged, i.e. Fi = —1.65, 
T2 = —2.05, 7 = 0.0), show how a different onset point on 
the MS leaves systematic deviations from the observed MS 
star counts just above and below the onset, caused by the 
(then wrong) MS lifetimes of the corresponding stellar mod- 
els. With hindsight, a mismatch in MS4, related to the use 
of an onset at 1.8 Mq, can be noted in our earlier synthetic 
samples (Schroder & Sedlmayr 2001), but the much smaller 
observed d < 50 pc sample used then as a reference did not 
provide a precise-enough test on the MS overshoot-onset. 

— Age of the thin disk population: We then made a 
comparison (again using the simple diffusion model) between 
synthetic samples computed for a different age of the "thin 
disk". These are listed in Table 2 as models t9ma X 10 and 
t9max8, where star formation starts 10 9 yrs earlier (see Bin- 
ney et al. 2000 who argue for an older solar neighbourhood) 
and later. The resulting synthetic samples are not very sen- 
sitive to the exact value adopted for age of the thin disk. 
This can be understood, since the contribution of the earli- 
est star formation to the present-day star counts is depleted 
by about a factor of 2. 

— Depletion description: To demonstrate the importance 
of matching the depletion factors, we also list a synthetic 
sample without any depletion (r<iif = oo), which models the 
solar neighbourhood (unrealistically) as a kind of "closed- 
box". Despite a significant reduction in the IMF of the 
less massive (older) stars to match present-day star counts 
(ri = —0.5 instead of -1.65), it is impossible to match 
the observed numbers of evolved stars. These are much 
lower than those predicted by stellar evolution time-scales - 
clearly, a depletion is required which progresses with age. 

See Fig. 2 for a plot of the depletion factor as a function 
of stellar age, compared to the reciprocal scale-height on an 
appropriately adapted scale (1 = 1/195 pc). Both curves 
coincide well for ages larger than 7 • 10 8 yrs, for which the 
dilution approach applies. 



in 



195pc/H_z 




2x10 9 4x10 9 6x10 9 8x10 9 

Age, in 1.E9 yrs 

Figure 2. Empirical depletion factors, averaged over star groups 
listed in Table 3, with approximate age ranges indicated. Open 
squares: depiction according to model (a), filled squares: model 
(b). Shown for reference are the depletion prescriptions for indi- 
vidual age: the exponential depletion used by model (a), and the 
dilution (<x -f/J 1 ) used by model (b) beyond the age required for 
thermalization and radial mixing (t* > 7 • 10 8 yrs). 

Finally, Table 2 lists several samples with the above- 
described approach of a depletion by dilution. The entry of 
Dil.(t* > 0.7 Gyr) represents the best match, obtained with 
Ti = -1.70 and T 2 = -2.10 for the IMF and without any 
change over time of the genuine SFR (7 = 0.0) in the sample 
volume. Table 3 and Fig. 2 show that the actual depletion in 
the various groups of stars of this model are indeed very sim- 
ilar to those derived from the simplistic but well-matching 
diffusion approximation. The final synthetic samples in Ta- 
ble 2 demonstrate the sensitivity of the synthetic samples to 
changes in the parameters for IMF and SFR. 

— Cut-off for z-values: We have also tested our anticipa- 
tion that it would matter to select for stars with small z 
when deriving the IMF from a given volume in the galactic 
plane. For this purpose, we modelled the stellar content of 
the whole spherical volume in d < 100 pc, i.e., with the larger 
2-values included. The best match (not listed here) then in- 
deed requires a much steeper IMF in its upper part, with 
r 2 = —2.35 (Ti = —1.75), caused by the already significant 
drop in massive stars with increasing \z\. 

— Best matching IMF and SFR: Table 4 summarises 
the properties of the IMF and SFR of the best matching 
synthetic samples. Despite their different approach to de- 
scribe the depletion of the galactic plane, both models (a) 
and (b) result in a quite steep IMF with Ti = —1.7 ± 0.1 
for 1.6M Q > M, >£s 1.1M Q and T 2 = -2.1 ± 0.1 for 
w 4M Q > M, > 1.6M Q . 

Furthermore, we find a present star formation rate 
SFR(0) in the solar neighbourhood of about 3.2(±0.1) ■ 
10~ 6 */yr in a volume of 1.5 • 10 6 pc 3 centred on the galac- 
tic plane, which is 2.1 stars per 1000 years and (kpc) 3 . This 
value of the SFR applies only to the stars with a mass larger 
than O.9M0. We cannot consider any less massive stars be- 
cause these are not luminous enough to appear in any of 
the specific HRD regions counted here, nor would they have 
evolved into the giant branch yet. 

In case there may have been a genuinely higher SFR 
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Table 3. Averaged properties (masses, ages, scale heights, depletion factors) with the simple diffusion approximation (a) and with stellar 
depletion by (mainly) dilution (b) for the synthetic stars of specific HRD regions (compare with Table 1). 





MSI 


MS2 


MS3 


MS4 


MS5 


MS6 


MS7 


KGC 


LGB 


CW 


< M > /M e 


4.0 


2.6 


2.0 


1.70 


1.53 


1.37 


1.24 


1.70 


1.30 


1.53 


<Age>/yrs 


1.9 10 8 


3.9 10 s 


5.9 10 8 


9.4 10 8 


1.43 10 9 


2.0 10 9 


2.8 10 9 


3.2 10 9 


5.1 10 9 


4.1 10 9 


< H z >/pc 


80 


135 


170 


210 


250 


280 


320 


350 


430 


390 


-Fdepl (a) 


0.97 


0.94 


0.91 


0.87 


0.80 


0.74 


0.68 


0.66 


0.48 


0.54 


fdcpl(b) 


0.96 


0.92 


0.87 


0.81 


0.75 


0.69 


0.63 


0.60 


0.46 


0.53 



Table 4. Slope of IMF (Scalo notation) listed for the two best 
matching models (simple diffusion approximation and depiction 
by diffusion & dilution, see Table 2), both with constant SFR, 
and two less good matches using a SFR slightly higher in the 
past. 



Depletion by 


SFR(0) 






r 2 


|<5|av 


T dif = 6.3 10 yrs 
Dil.(t* > 7 10 s ) 


3.03 10~ 6 
3.2 10~ 6 


0.00 
0.00 


-1.65 
-1.70 


-2.05 
-2.10 


0.30 
0.49 


Dil.(t. > 7 10 8 ) 
Dil.(t» > 7 10 8 ) 


3.0 10~ 6 
2.8 10~ 6 


0.03 
0.06 


-1.60 
-1.40 


-2.00 
-2.00 


0.72 
1.16 


in the past (9 billion years a£ 


;o) on 


a 1.3 


to 1.7x 


higher 



level (averaged radially and over a time-scale of a billion 
years), we tried models with 7 = 0.03 and 7 = 0.06, respec- 
tively. Such synthetic samples require a less steep IMF (see 
Table 4), but the over-all quality of the best match is sys- 
tematically less good (larger |<5 av |) than is achieved with the 
constant level of "thin disk" star formation. In particular, 
the count of the oldest group of stars (LGB) becomes signif- 
icantly too large (see last lines in Table 2) when 7 exceeds 
0.03. 

Finally, Fig. 3 shows the approximate star formation 
rate in the column (acc. to Equ. 15), as it appears from our 
observations in the volume in combination with our semi- 
empirical density scale-height as a function of age (Equ. 
10, 11). As already mentioned, the column integral SFR co i 
grows by a factor of 4 from its low present-day (local, inbe- 
tween spiral arms) value until it reaches an approximately 
constant, non-local disk-average of 0.82 * Gyr _1 pc~ 2 for 
all stars older than about 7 • 10 8 yrs (and M, > O.9M ). 
Of these, a large fraction must have been collected by the 
present-day solar neighbourhood as a result of radial mix- 
ing. We would like to note that a sampling range on the kpc 
scale (like the one for luminous stars of Scalo, 1986) could 
not show such a pronounced effect. Rather, a large sam- 
ple volume partly includes the neighbouring spiral arms and 
should show SFR C oi values much nearer to the disk-average 
in the first place. 



6 DISCUSSION AND OUTLOOK 

In our approach to assess the local galactic "thin disk" IMF 
of single field stars by means of synthetic stellar samples 
which match the observed, complete (volume-limited) sam- 
ple of solar neighbourhood stars, we have given special at- 
tention to three points vital to an unbiased IMF assessment: 
discrimination for stars near the galactic plane (\z\ < 25 pc), 
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SFR_col (M_* > 0.9 M_sun) 




8 8.5 9 9.5 10 

log(Age). in years 

Figure 3. Column-integrated SFR co i of the galactic disk over 
stellar age derived from our model (b), given in 10 — 9 stars per 
year and pc 2 , over log(agc). Note the apparent increase over the 
recent past (t* < 7 • 10 8 yrs) where radial mixing changes the ob- 
served SFR from its low, local (intcrarm) value to a much larger 
(and fairly constant) non-local average. For reference, the appar- 
ent SFR in the local volume is also shown (x200 pc). Its decline 
with age reflects the dilution of "thin disk" stars as they expand 
into larger scale-heights. 



a realistic (semi-empirical) quantification of the stellar de- 
pletion in the galactic plane, and use of the specific, age- 
dependent information buried in the account of the evolved 
stars in the stellar sample. 

Although we demonstrated that the simple approach of 
an age-independent diffusion time-scale, as used by us before 
(Schroder & Sedlmayr, 2001), does not yield too much dif- 
ferent an IMF from the more sophisticated, semi-empirical 
depletion by diffusion and dilution, the IMF derived here su- 
persedes our earlier work in several respects: (i) The present 
work is based on a much larger observed sample and yields 
significantly smaller error bars, (ii) The observed sample is 
a better representation of the stellar component in the im- 
mediate galactic plane, and (iii) the synthetic samples have 
now been created with a better choice of the onset of over- 
shooting on the MS. 

The slope of the field star IMF (i.e., F- L = -1.7, T 2 = 
—2.0) obtained by this study for medium to high masses 
is remarkably similar to the range of values considered by 
Scalo (1986, 1998), despite being based on a very different 
approach. In principle, the genuine IMF and the depletion- 
corrected SFR in a defined volume on the galactic plane, as 
determined here, are equivalent to their column-integrated 
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quantities studied by the classical approach, in how they 
depend on mass and age. In fact, in our column-integrated 
quantities, a depletion by simple geometrical dilution must 
cancel with the growing scale height. At least, this is true 
for stars older than 7 • 10 8 yrs. 

However, for stars younger than that, our column- 
integrated values are changing considerably, as demon- 
strated for the SFR by Fig. 3. This is the result of a pro- 
gressive radial mixing into the small, local and inter-spiral- 
arm volume used here. For the respective stars (mainly, 
M» > 2M Q ) our column- integrated IMF would be unphysi- 
cally steep. It does not compare to Scalo's (1986) IMF be- 
cause his sampling range is about an order of magnitude 
larger. Also, he uses different scale heights for that range 
of young stars. Still, upon closer inspection (Pagel, 1997, 
p. 207), a slightly steeper slope (F = —2.4) does actually 
appear in the respective data of Scalo (1986). This may be 
due to the same effect (radial mixing - hence, not genuine), 
but much reduced by his vastly larger sampling scale. In any 
case, the slope of the lower mass IMF remains unaffected. 

We see our approach of monte-carlo-created stellar 
samples, based on real evolution tracks and a good semi- 
empirical quantification of stellar depletion, compared one- 
to-one with an observed complete sample, as complemen- 
tary to recently developed methods of inverting colour- 
magnitude diagrams with a maximum entropy method (see 
Hernandez et al. 1999). While both methods make use of 
the full information content of the HR diagram, we find our 
approach is a very versatile tool and it may be more robust 
when studying the impact of more than one parameter. 

Still, the above-mentioned problems with radial mixing 
and a very limited sampling-range make it difficult to deduce 
an exact slope of the galactic disk IMF for the more massive 
stars. In fact, we are unable to adress the important question 
whether it may turn into a Salpeter dependence for the most 
massive stars (beyond about 6Mq) which are important for 
the UV radiation field and any interaction with the ISM. 
Rather, the size of our observed sample is more limited than 
we would wish for. This is also true when it comes to the 
question of any genuine change of the SFR, or even the IMF, 
over time (see discussion by Kroupa 2001, e.g. ). The results 
presented in the previous section provide, at best, support 
for an approximately constant SFR. 

But these limitations will disappear as the situation is 
going to change dramatically in about a decade, when a new 
dimension of data, in a comfortably large sampling range, 
will become available with the above-mentioned projects 
GAIA and DIVA. Well tested, detailed models for creating 
synthetic stellar populations will then be a very powerful 
means of assessing the star formation history in our galaxy, 
and the depletion by dilution which provides a sensitive test 
of our understanding of galactic dynamics. 
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